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1,  Introduction,. 

If  the  characteristic  function  of  an  n-dimensional  random  vector  x 
it  *  u 

has  the  form  e  -  -<J>(t'Et)  >  where  y:n*l  ,  E:nxn  ,  and  Z  >_  0  ,  we  say 
that  x  is  distributed  according  to  an  elliptically  contoured  distribu¬ 
tion  with  parameters  \ i,  Z  and  (j>  ,  and  we  write  x  ^  ECn(y,Ef<J))  (cf . 

Cam  ban  is,  Huang,  and  Simons  (1981)).  In  particular,  when  y  *  0  and 
Z  *  I  (the  identity  matrix),  ECn(0,In,<t>)  is  called  a  spherical  distri¬ 
bution  (cf.  Kelker  (1970)). 

The  class  of  the  elliptically  contoured  distributions  contains  such 
distributions  as  the  multivariate  normal  distribution,  the  multivariate 
t-distribution,  the  multivariate  Cauchy  distribution,  the  multivariate 
Laplace  distribution  and  the  multivariate  uniform  distribution  on  the 
sphere  in  Hn  or  in  the  sphere  in  Hn  . 

The  theory  of  elliptically  contoured  distributions  has  been  discussed 
by  Schoenberg  (1938),  Lord  (1954),  Kelker  (1970),  Das  Gupta  et  al.  (1972), 
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Devlin,  Gnanadeslkan  and  Kettenring  (1976),  Karlya  and  Eaton  (1977),  Muirhead 
(1980)  and  Cambanls,  Huang,  and  Simons  (1981).  The  following  basic  properties 
obtained  by  the  above  authors  are  needed  in  this  paper. 

(1)  Let  u  ^  denote  a  random  vector  which  is  uniformly  distributed 
on  the  unit  sphere  in  lRn  and  fi^dltlP)  denote  its  characteristic  func¬ 
tion.  If  ,  n  >  1,  is  the  class  of  all  functions  <J>:[0,°°)  K  such 
that  *(l|t|I2)  is  a  characteristic  function,  then  <f>  fc  if  and  only  if 

(1.1)  <J>(t)  -  f  0  (r2t)dF(r)  ,  t  >  0, 

Jo  a 

for  some  distribution  function  F  on  [O,00)  (cf.  Schoenberg  (1938)). 

(2)  x  ^  EC  (ll9£,40  with  rank  rk(Z)  «=  k  if  and  only  if 

(1.2)  x  -  y  +  FA'u(k) 

(k) 

where  R  ^  0  is  Independent  of  uv  '  ,  E  »  A* A  is  a  factorization  of  E 
(i.  e. ,  A  is  a  k  *  n  matrix  and  rk(A)  *  k)  and  the  distribution  func¬ 
tion  F(x)  of  K  is  related  to  <p  as  in  (1.1)  with  k  substituted  for 
n  •  Here,  x  *  y  denotes  that  the  random  vectors  x  and  y  are  identi¬ 
cally  distributed  (cf.  Cambanis,  Huang,  and  Simons  (1981)).  In  the  next 

,,dM 

section  we  discuss  the  operation  ■  in  detail. 

(3)  If  x  ^  Nn(y,Z)  ,  then  x  ^  ECn(y,£,$)  with  $(t)  •  exp(-t/2) 

2  2  2 

and  R  *v»  y  »  where  denotes  the  chi-squared  distribution  with  n 

degrees  of  freedom. 

If  A  is  an  n  x  n  symmetric  matrix,  what  is  the  distribution  of 
x'Ax ?  Kelker  (1970)  obtained  the  distribution  of  x'Ax/x'x  in  the  case 
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of  X  ^  EC  (0,1  ,40  and  the  distribution  of  xfAx  under  the  condition 

-  n  ~  -n  -  — - 

that  x  has  finite  fourth  moments  and  x  has  a  density.  Kariya  and 
Eaton  (1977)  gave  the  distributions  of  a'x/lixjl  and  x’Ax/I|x|P  ,  where 
||xll^  ■  x'x  .  There  is  a  rich  bibliography  on  the  distribution  of  quad¬ 
ratic  forms  in  normal  population.  Recently,  Anderson  and  Styan  (1980) 
reviewed  various  extensions  of  Cochran’s  Theorem  in  a  bibliographic  and 
historical  setting.  Our  purpose  In  this  paper  is  to  extend  Cochran’s 
Theorem  to  elliptically  contoured  distributions  in  various  aspects.  The 
main  results  are  in  Sections  4  and  5.  In  Section  3  we  list  some  basic 
properties  of  Dirichlet  distributions  that  are  needed  later.  Some  appli¬ 
cations  are  given  in  Section  6. 

Throughout  the  paper  N^y.Z)  denotes  the  n-dimensional  normal  distri- 

2 

bution  with  mean  y  and  covariance  matrix  E  ,  denotes  the  chi-squared 

distribution  with  k  degrees  of  freedom,  F(k,&)  denotes  the  F-distribution 

* 

with  k  and  l  degrees  of  freedom,  t^  denotes  the  t -distribution  with 
n  degrees  of  freedom,  I  denotes  the  n  x  n  identity  matrix,  denotes 
the  n  x  1  vector  with  elements  1  ,  rk(A)  denotes  the  rank  of  the  matrix 
A ,  and  A*  denotes  a  generalized  inverse. 

Md,f 

2.  The  Operation  =  . 

If  random  vectors  x  and  y  have  the  same  distribution,  we  denote 
that  fact  by  x  ■  y  . 

(1)  Assume  x,  y,  z  and  w  are  random  vectors,  x  and  z  are  lndepend 
ent,  y  and  w  are  independent,  and  x  -  y  .  Then  z  -  w  if  and  only  if 

(2.1)  x  +  z-  y  +  w. 
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In  particular,  if  z  is  a  nonrandom  vector,  then  x  ■  y  Implies  x  +  z 

d  A 
■  y  +  z  • 

(2)  Assume  that  x  *  y  and  f^(»)  ,  j  -  l,...,m,  are  Borel  func¬ 
tions,  then 


(2.2) 


fl(?> 


fMy)1 


f2(?>  d  f2(y) 


f(x) 

ID  ^ 


, .  d  , .  id. 

x  Ax  *  y  Ay  ,  x  x  ■  y  y 


xf  A^x 

fy'^y 

d 

xfA  x 

J 

I'hJL 
\  / 

For  Instance,  we  have 

(2.3) 

(2.4) 


(3)  The  following  fact  is  important  in  this  paper. 

Lemma  1.  Assume  that  x  and  y  are  n  *  1  random  vectors,  z  is  a 
random  variable  and  is  independent  of  x  and  y  ,  respectively.  If 

(2.5)  P(z  >  0)  -  1 

then  x  ^  y  if  and  only  if  zx  ^  zy  . 


Proof.  Firstly,  we  prove  that  Lemma  1  holds  if 

(2.6)  P(x,  >0)  -  P(y,  >  0)  -  1  ,  i-1 . n 


By  using  (2.2)  we  have 


Secondly,  we  consider  x  and  y  to  be  arbitrary  random  vectors.  It 
is  enough  to  prove  the  lemma  in  the  case  of  x  and  y  being  scalar  random 
variables  because 

x  -  y  °  afx  -  afy  >  Va  fclR  **  zafx  -  zafy  ,  Va  €lR 

a.  x#  -  ^  <v  n  ■v  *v  —  — ■  n 

**  a'(zx)  ■  a'(zy)  .  Va  fc]Rn .  *»  zx  ■  zy 

If  we  have  proved  the  lemma  in  the  above  case.  Let 

x  if  x  >  0  0  If  x>0 

f+(x)  -  <  ,  f_(x)  -  - 

0  if  x  <  0  [-x  If  x  <  0 

and  x+  “  f+(x)  ,  x”  ■  f_(x)  ;  y+  *  f+(y)  and  y  ■  f_(y)  .  It  is  easy  to 
verify  that 

x  ^  y  •  x+  ■  y+ ,  x”  ■  y”  ,  and  p(x*0)  ■  p(y*0)  *  x+  ■  y+  and  x  ■  y  . 

If  we  can  prove  that  zx+  ^  zy+  »  x+  ^  y+  and  zx"  ^  zy"  «*  x"  4  y"  ,  then 
the  essertlon  follows  from 
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zx  ■  zy  •*  (zx)+  •  (zy)+  and  (zx)~  *  (zy)”  «* 

+  d  +  ,  -d  -  _  +  d  +  .  -  d  -  „  d 

zx  •  zy  and  zx  •  zy  **  x  ■  y  and  x  *  y  »  x  •  y . 

.  .  +  d  +  _  +  d  +  A  *f  d  +  , 

Now  we  prove  that  zx  ■  zy  *  x  *  y  .  Assume  zx  *  zy  ,  then 

p(x>0)  ■  p(x+>0)  ■  p(zx+>0)  ■  p(zy+>0)  ■  p(y+>0)  “  pCy^O) 

and 

+  + 

p(x>0)E(e*tZX|x>0)  “  E(eitZX  )  -  p(x<0)  *  E(eitzy  )  -  p(y<0) 

-  E(eitzy|y>0)  p(y>0). 

From  the  first  part  of  the  proof  we  have  E(e*tx|x>0)  *  E(eity|y>o)  , 
i.e. ,  x+  ■  y+  .  The  part  follows  by  the  same  technique.  Similarly 
we  have  zx  ■  zy  **  x”  *  y”  .  Q.E.D. 


3.  The  Dirichlet  Distribution. 

If  y  »  (y^, . . .  .y^)  '  is  a  random  vector  with  E™  y^  *  1  and 
(y^» • • • ,ym  ^) *  has  the  density  with  >  0  ,  i  *  1, . . . ,m  , 


(3.1) 


pm^tl*“*,tm-l^ 


m 

b 


ra-l  a 
n  t 


i-l 


nr  (a.) 
i  1 


m-l 

>  It. 
1 


a  i 

m— 1 


if  t^  ^  0  ,  ,m-l  , 


1 


■  0  ,  otherwise. 

We  say  that  y  is  distributed  according  to  a  Dirichlet  distribution  and 
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I  t— J«H 


When  m -  2  D2(d1;d2) 


denote  it  as  (yx . y^)  *  . Vl;CXm)  * 

reduces  to  the  Beta  distribution,  B(a^,a2)  . 

It  is  a  well-known  fact  that  if  x  ^  N  (0,1  )  and  is  partitioned 

—  n  —  ~n 

into  m  parts  X(x) > *  *  * *x(m)  with  ai»***>aa  components  of  x,  respec¬ 
tively,  i.e. , 


(3.2) 


?(1) 


(m) 


then  the  joint  density  of 


-(!)?(!) 


.I!2  . 


X  *  X  ^ 

~(m-l)~(m-l) 


is  D 
m 


a  ,  a 
m-1 .  m 

»  2  *  2 


In  particular,  if  a  *...  “0^*1,  Qt^+^  ■  n  -  k  ,  m  «  k+1 ,  then 

has  Djt+i  •  •  •  »y  »  as  its  density.  It  can  be  veri- 


2  2 
i  xk 

,2,‘"’l|x||2. 


[h\\ 

fled  that  the  density  of 


«*1 1  »*k 


~  ) 


n-k 


(3.3) 


r(|)  2k  f 

r(5— )iri5k 


-  1 


1  -l  -l 
1 


is 


if  x 


k  2 

>0,  i“  1,.  • .  ,k  ;  ^X-<1, 

1  “  1  x 


and  the  density  of 


1  k 

[W . WJ 


is 


(3.4)  8q(x^,.o,x^)  — 


r(f) 


rc 


"iV71 


n-k 


-  1 


i  - 1  *i 

i  * 


,  if  i  xi  <  1 

i  1 


If  x*  E  C  <0,1  W  ,  then  (cf.  (1.2)) 
*»  n  -  •**n 
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where  R  end  uv  ' 


are  Independent.  From  (2.3) 


(3.6)  j|*Jp  -  x'x  »  R2u^  “  R2  or  ||x||  ^  R  . 

If  p(x*0)  ■  0  (or  p(R“0)»0) ,  we  have  (cf.  Cambanis,  Huang,  and  Simons, 

•» 

(1981),  Lemma  1) 

(3.7)  x/{[x |{  -  u(n>  , 


and  it  is  independent  of  ||x||  ^  R.  Note  that  (3.7)  holds  for  all  of  the 
class  of  spherical  distributions  with  p(x“0)  *0.  Thus  we  can  assume 
x  ^  N  (0,1  )  to  obtain  the  distribution  of  ratios.  We  will  use  the  property 

n  -v  -n 

many  times.  Denote  u^  -  (u^, . . .  ,uQ)  * .  We  Immediately  obtain  that 
(i)  the  density  of  u.  , . .  .  ,u  ,  1  £  i.  <  i9  <  •  •  •  <  i.  £  n  is 

h  *k  1  1  * 

.  ,Xj^)  ,  (cf.  (3.4)); 

(ii)  the  density  of  |uA  | , . . . ,  ju^  | ,  1  <_  i^  <  i2  <  •••  <  ifc  <_  n,  is  (3.3) 
f  \  ^  k 

(Hi)  if  u',n‘'  is  partitioned  into  m  parts  as  the  same  as  (3.2),  then 
the  density  of  (||u(1)!|2 . llH(*-i)ll>  is  Vi  al . T  Vl;  T  am)# 

tn 

(Iv)  Let  y  ■  E[Hu.  ]t  then  y  *  0  if  some  r . 

rl*  *  *  *  ,rm  1  1  rl*  *  *  *  ,rm  1 

are  odd.  If  all  of  r-,..«,r  are  even,  y  can  be  obtained  as 

i  m  rl,***,m 

moments  of  the  Birichlet  distribution  (cf.  Johnson  and  Kota  (1972)),  i.e. 


m  i  [riJ 

n  ih  j 


(3.8) 


m 

II  r,] 

tfi1  J 


where  x^  m  x(x+l)  •••  (x+k-1) .  in  particular,  Eu^  -  0,  Bu^V^  ■  —  I 


There  Is  a  dose  relationship  between  the  chi-squared  distribution  and 
Dirlchlet  distribution  from  the  above  discussion.  The  following  are  further 
results. 


Lemna  2.  Assume  that  (z^  . . .  ,zm>  ^  aml  j  ot^)  with 

>  0,  i  ■  l,...,m+l,  and  *  n  and  yo,yi*  •  •  •  ,ymji  are  distri¬ 
buted  as  chi-squared  distributions  with  degrees  of  freedom  n,a. , . . . ,a :  , 

1  nH-1 


respectively;  then 


(3.9) 


h  z,  ,...,£  z  (t1»‘**»tm) 
n  X  n  n 


1  V,(v 
%>oi  v 


r(f)  m  rc^+it^) 

r(f+i(t1+...+tB))  j-i  r(^L} 


where  $(•)  denotes  the  characteristic  function  and  i  *  v^T 


Proof.  Let  x  ^  NR(0,In)  and  x(i) * • • • »x(m+i)  k*ve  similar  sense 
as  (3.2).  Then  fcojll*  *  y j ,  j  *  and  IWI2  -  y^.  We  have 


ll;(1)ll2l  fll;(l)ll2/l;l2] 


U;<l0ii2J 


[‘Jji 


2  ill  if2^ 


i  ‘-Mi2  in  * 


l‘JI«(.,ll2/L«ll21 
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are 


The  first  part  of  (3.9)  follows  from  the  fact  that  ||x||  and  x/j|x|l 

independent  and  (l;(1)ll2/«;ll2 . Il;(»)ll2/ll;ll2>  ■  (Zj^ . Z^)  .  Note 

(cf.  Press  (1969)) 


itf  y  2ilT(-f+it) 

E(e  J) - L -  ; 


r£> 


the  second  part  of  of  (3.9)  follows.  Q.E.D. 


In  particular,  if  z  ^  B(-|»  ^2^),  yQ  ^  and  ^  X^»  then 


\yl(t)  r(f)r(£+it)  B(f+it,  ^) 

(3*10)  ♦*.  z(t>  “  ?T7uT  -  r(|)r(f+it)  * 


n 


Vo 


2’  2 


where  B(a,b)  is  Beta  function. 


Lemma  3.  Assume  that  x^N  (0,1),  (z  ,  ...,z  .)  'v  D  (L..  ~;i), 

~  n  ~  ~n  x  n-x  n  z  z  z 

2 

zq  ”  1-1^  z^,  y^,...  ,yn  are  i.i.d.,  y^x^  1  “  1 . n,  and 

A^,...,Ajn  are  n  x  n  symmetric  matrices.  Let  gi(t^, . . .  ,t^)  ,  i  -  1 . . 

be  linear  functions  of  t, ,...,t  ,  then 

1  n 


(3.11) 


Vh;/n?n2' 

1 

Vsi . *„>’ 

* 

d 

• 

• 

• 

• 

• 

gm(z1,...,zn)j 

ft  y 

x'A^x1 

. V' 

• 

d 

• 

• 

m 

| 

i 

x'A  xl 

8m(yl,*“,yn)> 

Proof.  We  firstly  prove  the  implication  towards  the  left.  If  the  right 
side  of  (3.11)  holds  then 
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If 


VjjX/llxf 

• 

x '  Aj^x 

• 

f 

d 

• 

• 

• 

• 

m 

*  *  *  »xn) 


2 

‘l* 


^(x, , . . . , 


Xn>i 


INI2 


8i<: 


2 

X1 


'W  MT 


. IUII2’ 


and  from  Lemma  1 


V^x/llxf- 

• 

d 

■ 

2  2  i 

X1  x 

8l<iuF . wf’ 

d 

m 

8l(zl,,,,,zn) 

• 

• 

• 

x'Ax/l|xl|* 

2  2 

X1  Xn 

Hxif . n*r 

• 

V*1 . V' 

2 

Here  we  use  the  fact  again  that  ||x||  and  x/1|x||  are  independent.  The  proof 
in  the  other  direction  is  similar.  Q.B.D. 

By  Lemma  3,  we  can  change  the  theory  of  distributions  of  quadratic  forms 
from  chi-square  distributions  to  Dirichlet  distributions. 

Corollary  1.  Assume  that  x  *v»  Nn(0,I)  and  A  is  an  n  x  n  symmetric 
matrix,  then  x'Ax/|lx|p  'V-  B(y,  if  and  only  if  A^  ■  A  and  rk(A)  ■  k. 

2 

Proof .  By  Cochran's  Theorem,  A  ■  A  and  rk(A)  ■  k  if  and  only  if 
x'Ax  xl  or  x'Ax  ^  hence  if  and  only  if  x'Ax/||x|f  ■  ^  B(y, 

from  Lenana  3  with  m  ■  1  and  •  •  •  ttn>  *  ti+***  +  tk’  Q*E,D‘ 
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Corollary  2.  Assume  that  x  ^  N_(O.I_),  A  and  B  are  n  x  n  symmetric 

1  —  ■■  1  *  -  q  -v  ~n  «  •« 

matrices,  then  z. . z  .  has  the  Dirlchlet  distribution  D  (4- , . . . ,  ^ ; -I)  and 

x  u*i  n  1  2  a 

Ej  «4  ■  1 .  then 


(3.12) 


x'Ax 

'0 


k+m 

l 


k+l 


V  z 

j  j 


k+m  £  n  , 


If  and  only  if  AB  •  0, 


where  A  ^  and 


are  nonzero  real  numbers. 


Proof .  By  Craig's  Theorem  (cf.  Anderson  and  Styan  (1980))  AB  *  0  if 
and  only  if  x'Ax  and  x'Bx  are  independent  or  if  and  only  if 


(x’Ax,  x'Bx) 


II 


Vi’ 


|  Vj*j 


where  X, ,...fX,  and  v  ,  are  the  nonzero  eigenvalues  of  A  and  B, 

lk  1  m  ~ 

respectively.  Thus  the  corollary  follows  by  Lemma  3.  Q.E.D. 


Corollary  3.  Let  x,  A  and  B  be  defined  as  in  Corollary  2,  then 
(x* Ax/|[x|f ,  x’Bx/Kxli2)  *  Dj(y,  ,  with  k  >  0,  m  >  0,  k+m  <_  n,  if 

and  only  if  AB*  0,  A2  -  A,  B2  -  B,  rk(A)  -  k  and  rk(B)  -  m. 

k  k+m 

Proof.  The  corollary  follows  from  Corollary  2  and  (E^  z^,  E^+^  z^)  ^ 

_  .k  m.  n-kmx  _ 

DjIy*  2*  2  *  Q.E.D. 

4.  Distributions  of  Quadratic  Forms  and  Cochran's  Theorem. 

In  this  section  we  want  to  extend  distributions  of  quadratic  forms  and 
Cochran's  Theorem  to  the  case  of  elllptlcally  contoured  distributions  with 
mean  zero.  Firstly,  we  need  to  generalize  the  Dirlchlet  distribution  and 
the  Beta  distribution. 
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If  a  random  vector  (z, )'  satisfies  (z.,...,z ;  )  -  R  (u. 

X  B  X  m  XI 

where  R  ^  F(x),  R  is  independent  of  (u^, . . . *un_^) »  2®  ■  1,  and 

(ulf. ..,Ub1)  -V-  then  we  write  (Zj, . . ^ 

G  (a_  , . . .  ,o  -  }<x  ;<^)  and  (z.  t*>«,z  )  ^  G  (ct. ,  •  • .  ,oi  .  ,ci where  <f> 
a  i  n-i  m  x  m  m  x  m-x  m 

is  related  to  F(x)  as  in  (1.1)  with  n  ■  2(a,  +  •••+<*  ). 

x  in 

It  is  easy  to  show  that  the  density  of  G  (a. . a  .  *a  ;<t>)  is 

m  x  m-x  in 


(4.1) 


r(f) 


n  r(o  ) 

i  1 


B-l 

n 

i 


V1 


-(n-2) 


/m-1 

/  J«i 


m-1 

-  I « 

i 


a  -1 
m 


dF(r)  ,  if 


Further,  if  R  has  a  density  f(r),  then  the  joint  density  of  u^,...,u 
and  R  is  (cf.  (3.1)) 


r<!>  "-1  v1 ' 

-  n  u 

i=l 


n  r(a.) 

i-l 


n-1 


a  -1 
m 


f  (r) 


Consider  the  following  transformation: 


z±  —  r  u^,  i  ■  1 , • • • »m-l  , 


\  m-1  ) 

- 


then 


zJt  ,  i-l,. 


. ,m-l  , 


IM 

l  *i 


1/2 


Thus  the  Jacobian  of  the  transformation  is 
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n/2 


(4.3) 


n  r(a.) 

1 


m  o.-l 

n  «*  g 

i-l  1 


m 

h 


In  Section  6  ve  will  give  some  examples  of  G  (a, , . . .  ,a  ,  ;a  ;<t>)  .  Anderson 

m  X  m-1  n 

and  Pang  (1982)  give  some  further  applications  of  G  (a,,..., a  ,  ;a  ;<J>) . 

m  x  m-i  n 


Theorem  1.  Suppose  that  x  ^  EC  (0fI  »♦),  p(x«0)  ■  0  and  A  is  an 
nxn  symmetric  matrix,  then  x'Ax  ^  G^C^*  y^;  4>)  If  and  only  if  A^  -  A 


and  rk(A)  “  k. 


Proof.  If  A*  •  A  and  rk(A)  ■  k,  then  there  exists  an  orthogonal 


matrix  T  such  that  PAT1 


ik° 
0  0 


Let  y  *  T'x;  it  can  be  verified  that 


r.k  2 


that  y  •v  EC  (0,1  ,$)  and  x'Ax*  I.  y,,  where  y*  (y. , . . .  ,y  )  ' .  Thus 
n  ~  i  l  %  i  n 

k  „ 

9  ..9 

(4.4) 


t  ,  o  i  y  ^ 

x'Ax  -  I  y2  -  ||y||2  •  J-4 
~  ~  1  1  II  yli 


As  Is  known,  Hyl)2  and  E*  y2,t{\y\^  are  Independent  and  the  distribution  of 
9  1  1  ^ 

if-  y^/IMI2  does  not  depend  on  what  the  distribution  of  y  is  In  the  class 

1  1  'W 

of  EC  (0,1  ,$);  therefore  we  can  assume  y  N  (0,1  ).  It  implies  that 
n  %  — n  ~  n  -  — n 

E*  %  y^)*  As  l|y||  *  R  ^  F(x),  we  have  completed  the  proof 

of  the  "if”  part. 

If  x'Ax  'v*  G2(y»  y*;4>),  write 

,'ta  -  Nf  . 

-  “  ’  Ml2 


where  R  ^  P(x),  z  *  B(y,  y^)  and  they  are  independent.  By  Lemma  1  we  have 


15 


x'Ax 

IwF 


i 


z 


From  Corollary  1  of  Lemma  3  the  assertion  follows. 


Q.E.D. 


Corollary .  Assume  that  x  ^  ECn(0,E,4>)  with  I  >  0  and  p(x«0)  ■  0, 

It  pj— 

and  A  is  an  n*n  symmetric  matrix;  then  x’Ax  'v*  G2 (^S  —£~ ;  <J>)  if  and 
only  if  AEA  *  A  and  rk(A)  *  k. 

Proof.  As  E  >  0,  there  exists  E*4  >  0  such  that  E^E*4  ■  E.  Let 
y  ■  E  J*x,  then  y  'v  ECn(0,In,<|>) .  By  Theorem  1 

x’Ax  -  y' (E^AZ^y  'V-  G2^ 

if  and  only  if  (E^AZ*4)2  «  Z^AE*4  and  rkCE^AE*4)  -  k,  i.e.  AEA  -  A  and 
rk(A)  -  k.  Q.E.D. 

2 

It  is  well-known  that  A  ■  A  if  and  only  if  rk(A)  4-  rk(I~A)  *  n. 
Hence  AEA  “A  if  and  only  if  rk(A)  +  rk(E-A)  *  n.  The  corollary  shows 
us  that  E  must  be  a  generalized  Inverse  of  A.  Kelker  (1970)  gives  a 
result  similar  to  the  corollary,  but  he  assumes  finite  fourth  moments  for 
the  components  of  x. 

If  x  ^  N  (0,1  ),  A'  ■  A,  A^  *  A  and  rk(A)  ■  k,  then  x'Ax  ^  Xu  by 
Cochran's  Theorem,  i.e.  g2(j;  Sjp;  <|>)  -  x£  if  4>(t)  -  exp(-t/2).  The 
Inverse  proposition  also  holds. 

Theorem  2.  The  distribution  ^5^*  ♦)  is  x£  if  and  only  if 

♦(t)  -  exp(-t/2) . 


16 


Proof.  We  only  prove  the  "only  if"  part.  If  z  'v-  G2(-^l  <fr)  is 

x£,  by  the  definition  of  We  **ve  z  ^  R  ‘w,  where 

R  <v  F(x),  w  *v  B(j,  and  R2  is  Independent  of  w.  Let  yA  %  x£ 

2 

and  y^  *v*  x^»  from  (3.10)  and  the  supposition  we  find 


h  y  <C)  “  h  z(t>  “  \  R2(t)  •  V(t>  “  **  R 
n  1  n  *  R  n  *  R 

n  n 


♦t  y  <C> 
(t)  ~n  l 
2  v  (t>  * 

Vo 


which  implies 


i .  e . 


♦<t) 


exp(-t/2) . 


Q.E.D. 


By  a  method  similar  to  those  used  in  proving  Theorem  1  and  Theorem  2, 
we  obtain  the  following  theorems  and  corollaries. 


Theorem  3.  Suppose  that  x  ^  EC_(0,I_,<J>)  ,  p(x-0)  «  0,  A  and  B  are 

-  "  n  -  ~n  ~  ~  ~  ~ 

k  Z  n-k-Z 

n  x  n  symmetric  matrices,  then  (x'Ax,  x’Bx)  **•  G^O^,  -j*  — 2 — *  ^  **  and 

only  if  AB  -  0,  A2  «  A,  B2  -  B,  rk(A)  -  k  and  rk(B)  -  i. 

Corollary  1.  Suppose  that  x  ^  ECn(0,E,<|>)  with  E  >  0,  p(x-0)  -  0, 

k  Z  n-k-Jl 

A  and  B  are  n  *  n  symmetric  matrices,  then  (x'Ax.  x'Bx)  'v  G3(j,  y;  — 2 — » 
if  and  only  if  ALA  -  A,  BEB  -  B,  rk(A)  -  k,  rk(B)  -  l  and  AEB  -  0. 

and  B  are 
AB  •  0,  then 


Proof.  As  AB  -  0,  there  exists  an  orthogonal  matrix  T  such  that 
r’AT  -  diag(l,...,l,0,...,0)  and  r’Br  -  dlag(0, . . . ,0,1, . . . ,1.0, . . . ,0) . 
k  k  l 


Corollary  2.  Assume  that  x  ***  ECn(0,In,<|>) ,  p(x»0)  “0,  A 
projection  matrices  with  rk(A)  ■  k  and  rk(B)  *  1  satisfying 


x'Ax 


k  x'Bx 


*  F(k,L)  . 
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r 


Let  y  ■  Tx, 
thus 


Ch<n  y  *  ECn(0,In.«),  J-A*  .  y\  .nd  x'Bj  -  y\. 


i 

k  x'Bx 


£  1 


k  k+£ 


k+1 


k 

I 

£ _ 1_ 

k  kf£ 

I 

k+1 


y^iiyi 


^/llyll2 


k 

I 

d  £  1 
“  k  k+£ 

I 

k+1 


*  ^11 


Si/fl*U2 


£  1 


k  k+£ 


k+1 


*  F(k,£)  , 


where  z  *>»  N  (0,1  ).  Q.E.D. 

^  n  -  -n 

k  2  lr4.0  ? 

Kelker  (1970)  treated  the  distribution  of  Ej  Xi/E^+i  Xi  by  a  different 
method. 


Theorem  4.  The  joint  distribution  G^(yt  y;  ,n  2"  »  M  1®  the  product 
of  the  distributions  of  x£  and  if  and  only  If  <J>(t)  *  exp(-t/2)  . 


Theorem  5.  Suppose  that  x  ^  ECn(0,In»<fr) ,  p(x**0)  -  0,  A  and  B  are 
n  x  n  symmetric  matrices  and  (z^,...,z  .)  *\»  D  (y, ... ,4;  y)  Is  Independent 

of  R.  Then 


x'Ax 

x'Bx 


k 

pi*i 

k+£ 

l  Vi 

k+1  2  3 


with  real  numbers  X.  and  v.  if  and  only  If  AB  -  0  and  X .  's  are  the  non** 

X  j  —  X 

zero  eigenvalues  of  A,  v  's  are  the  nonzero  eigenvalues  of  B. 

~  J  ~ 

The  above  theorems  can  be  generalized  to  the  case  of  several  quadratic 
forms. 
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5.  Tripotent  Matrices. 

3 

A  square  matrix  A  is  said  to  be  tripotent  whenever  A  -  A.  Anderson 

and  Styan  (1980)  extended  Cochran's  Theorem  to  tripotent  matrices.  If 
3 

A  -  A,  the  eigenvalues  of  A  are  1,  -1  and  0.  Let  p  and  q  denote 

the  number  of  eigenvalues  equal  to  1  and  -1,  respectively.  If 

x  'v  N  (0,1  ),  then  x'Ax  ^  y.-y,,  where  y.,  and  y  are  Independent, 

2  2 

y^  'V  Xp  and  y2<vXq*  The  distribution  of  y^-yj  has  been  given  by  Press 
(1969)  and  Robinson  (1965)  and  we  denote  it  by  H(^,  -^) .  Similarly,  if 
(v^.v^)  D3^»  2*  — 2~^)  we  denote  the  distribution  of  z  ■  v^-v^  by 

H(f,  f;  S=P>. 


Lemma  4.  Assume  that  *  ^  N_(0,I_)  and  A  is  an  n  x  n  symmetric  matrix. 
-  n  -  ^.n  — 

Then  x’Ax  H(^-,  if  and  only  if  A^  ■  A  with  p  l's  and 

q  -l's  as  its  nonzero  eigenvalues. 

Proof.  Apply  Lemma  3  with  g^(t^, . . . ,tn)  ■  t^  -  tj  and  m  ■  1 

which  completes  the  proof. 

Now  we  generalize  the  distribution  H(^,  P~P~9)  to  the  case  of  elllp- 

tically  contoured  distributions.  If  z  ■*  vi“v2  ^  "f*  D~2'  ^  ®  ^ 

2 

and  is  independent  of  z,  we  denote  the  distribution  of  R  z  by 
H(|,  -P~P~9;  ^  where  $  is  related  to  F(x)  as  in  (1.1).  By  the  same 

technique  used  in  the  proof  of  Theorem  1  and  Lemma  4  we  obtain  the  following 
theorem. 

Theorem  6.  Suppose  that  x  ^  EC_  (0, 1_  ,<t>) ,  p(x*0)  -  0  and  A  is  an 

1  1  1  iv  Q  v  vll  v  v 

n  x  n  symmetric  matrix,  then  x'Ax  ^  H(^,  n~P~^;  $)  if  and  only  if 

A*  A  with  p  l's  and  q  -l's  as  Its  eigenvalues. 

«*• 
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Corollary*  Suppose  that  x  ^  EC^tO^S,^)  with  Z  >  0,  p(x-O)  -  0, 
and  A  is  an  n  *  n  symmetric  matrix,  then  x'Ax  ^  H(^-,  Q~P~<i .  ^ 
if  and  only  if 

(5.1)  rk(A)  -  rk(A-AZA)  +  rk(A+AlA) 


and  AE  has  p  l's  and  q  -l's  as  its  eigenvalues. 

Proof .  As  £  >  0,  there  exists  E*5  >  0  such  that  E**!**  ■  E.  Let 

x  •  eSt,  then  y  *v>  EC  (0,1  ,<f>)  and  x'Ax  *  y' (E^AE^y.  From  Theorem  6 

—  —  —  **  n  -  -n  ~  ^ 

x'Ax  ^  H(£,  —  f- 0)  if  and  only  if 

**  4  4  4 

(i)  Z^AZ  has  p  lrs  and  q  -l’s  as  its  eigenvalues,  or  equivalently 
AE  has  p  l's  and  q  -l's  as  its  eigenvalues* 

(ii)  (E^AZ*5)^  “  (E^AE**)  .  Anderson  and  Styan  (1980)  point  out  that 

o  2  2 

B  -  B  if  and  only  if  rk(B)  -  rk(j$4B  )  +  rk(B-B  )  for  any  square  matrix 
B.  Thus  (E^AI*5) 3  -  E^AE*5  if  and  only  if 

rk(A)  •  rkCE^S  -  rk  ( E^AZ^+E^AZ  AE**)  +  rk  (E^AE^-E^AEAZ*4) 

-  rk(A+A!A)  +  rk(A-AEA)  .  Q.E.D. 

6.  Applications. 

In  this  section  some  applications  are  given.  In  the  first  part  we 
consider  the  theory  on  linear  regression  with  the  error  being  distributed 
according  to  an  elllptlcally  contoured  distribution.  We  will  find  the 
distribution  of  the  sample  variance,  the  distribution  of  the  ratio  of  the 
sample  mean  to  the  sample  standard  deviation  in  the  second  part.  In  addition 


we  will  give  the  distributions  of  quadratic  forms  for  the  multivariate 
t-distributlon,  the  multivariate  uniform  distribution,  etc. 


We  consider  the  following  model: 


?n*l  "  ?n*p~pxl  +  ?nxl*  -  P  <  n 

(6.1)  <  E(e)  -  0,  E(ee')  -2>  0 

e  ^  ECn(0,E,^)  . 

There  exists  a  pxp  matrix  A  such  that  A'A  -  I  and 

(6.2)  e  ^  »A’u(n)  . 

Hence 

(6.3)  y  =  Xg  +  RA'u(n*  , 

i.e.,  y  ^  EC  (Xg,E,<f>)  . 

—  n  ~~  - 

Minimizing  e'e  ■  (y-Xg) ' (y-Xg)  with  respect  to  g  gives  the  least 
squares  estimator 

(6.4)  g  -  (X’X)X'y  . 

In  order  to  find  the  distribution  of  g  we  need  the  following  leoana. 


Leona  5.  Assume  y 'V  EC  (u.I  .d>)  and  A  is  n  x  o  matrix  with 

■  ~  n  -w  — n  i 

rk(A)  ■  p  <  n;  then  X  3  A'y  ^  EC  (A'y.A'A,^) »  where  <p  e  $  corres- 

~  ~  ~~  p„~~~  p 

.  * 
ponds  to  R  and 
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(6.5) 


R*  ^  R  •  b 


t 


2  11 

where  b  0,  b  *x»  B(-jP  ,  -|(n-p)) ,  end  b  is  independent  of  R. 

Proof.  Consider  the  singular  value  decomposition  of  A,  i.e.,  there 

exist  orthogonal  matrices  I*  and  A  _  such  that 

~PXP  ~nxn 

A'  -  r(Dr  0)A'  , 

where  ■  diag(X^ . X^)  and  X^,...,Xp  are  the  roots  of  |a'A-Xi|  ■  0. 

From  the  assumption 

y  ^  y  +  Ru^  , 

thus 

x  i  A'p  +  RA'u(n) 

-  A’y  +  RT(D^  0)A'u*n) 

^  A'y  +  RT(DX  0)u(n)  . 

because  u^n^  ■  A'u^n\  Let  z  'x*  N  (0,1  )  and 
-  ^  ^  ~  n  -  ~n 


r  ' 

) 

(n) 

u  ■ 

5(i) 

,  2  * 

~(D 

— 

?<*), 

.  -<2\ 

where  u/1%  and  z/1%  are  p  x  1  vectors.  Based  on  the  relationship  between 
-Cl}  -\l} 

u<n>  and  z  we  have 

X  i  A,y  +  Rrpxu(1)  i  A'y  +  I^rpx*(1)/|Jz|| 

•  ■<!.>  IWI  *  H>*  ( 1)  'll  *ci)  ® 

^  A'y+(Rb)(rD.)u<p)  . 

—  —  <v^A  "• 
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follows.  Q.E.D. 


Sine*  A’ A  -  (rD.)(rD1)',  the  Leona 

By  Lama  5  and  (6.4)  we  have 

(6.6)  §  “  + 

*  0  +  R(X'X)“lX'A,u^  , 

2 

where  4>  is  defined  by  Lemma  5.  In  particular,  if  E  «  a  I  ,  then 

~P 

0  *  ECp(0,a2(X,X)"1,<J>)  . 

Denote 

(6.7)  s  -  (y-X0) ’ (y-X0) 

-  y,[l-X(X’X)~1X,]y 

«S»  ~  •'S'  ~  *w  %  V 

»  e '  I  I-X(X'X)  "^X* )  e  . 

If  e  EC  (0,a2I  ,(j>),  p(e»0)  ■  0  and  rkx  ■  p,  then  s  *v-  o2G,(^E-;  •£»  ♦) 
-v.  n  -  -v  a  ^  a 

by  Theorem  1.  We  summarize  the  above  results  in  the  following  theorem. 

2 

Theorem  7.  Assume  e  ^  EC  (0,0  I  ,i),  p(e«0)  *  0  and  rkx  *  p  <  n, 

■  ■  «  n  ~  <vti  *v  v 

then 

0  *  ECp(0,O2(X,X)_1,().) 

and 

s  ^  a2G2(~E;  <f>)  . 


r 


Now  we  consider  the  linear  hypothesis 
(6.8)  H:  HB  «  c  , 

where 

H:  qxp,  rk(H)  *  q  <  p  . 

Under  the  condition  of  H$  *  c  minimizing  e'e  with  respect  to  8 
gives  the  least  squares  estimator 


(6.9) 


eH  -  §  -  (x'x)-1h'(h(x'x)-1(h3-c)  , 


where  8  is  given  by  (6.4),  and  the  corresponding  (y-XB^) ' (y-XB^) 
becomes 

(6.10)  sR  =  (y-X§H) 1  (y-x§H)  - 

-  s  +  (H$-c) ,(H(X’X)~1H,)~1(H6-c)  , 

where  s  is  defined  by  (6.7), 

The  statistics  testing  the  hypothesis  (6.8)  is 


(6.11) 


Vs 


When  the  hypothesis  is  true  we  have  c  *  HB  and 

(6.12)  (H8-c)  »  H(§-B)  -  Wi(X,X)"1X,A,u(n)  -  H(X,X)"1X,e  . 


■V  "v  *v 


Thus 


(6.13)  Sfl-s  -  «'t?(?'?)~1H,(H(X,X)"1H,)“1H(X,X)“1X,]e 


=  e'D^e  (say) 
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It  is  easy  to  see  that  pjj  -  pH,  rkpH  -  q,  and  pH(I-X(X’X)-1X')  -  0.  From 

(6.13)  and  (6.7)  X  can  be  expressed  as 

e%  e 

(6.14)  X - - - -  , 

e1  [I-X(X'X)  Xr]e 

and 

X  ^  F(q,n-p)  . 

Q 

(cf.  Corollary  2  of  Theorem  3). 

2 

Theorem  8.  Assume  e  ^  ECn(0,a  1^,$),  p(e~0)  *  0  and  rkx  =*  p,  when 
the  hypothesis  (6.8)  is  true  the  distribution  of  (n-p)A/q  is  F(q,n-q). 


6.2.  Some  examples. 


Example  1.  If  x  ^  EC^(0,I  ,<(>)  and  p(x«*0)  -  0,  let 
-  n  ~  -n  ~  — 


(6.15)  x  •  ~  l  x.  and  s 2  *  £  (x  -x)2  . 

1  1 

2  1 

We  can  express  s  *  xf(l  -M)x  with  M  *  —  £  e*.  It  is  easy  to  check  that 

~  -  —  ~  n  -n~n 

(I-M)2  *  I-M  and  rk(I-M)  *  n-1.  By  Theorems  1  and  2  s2  ^  G20^r^;  4>) 

and  s2  •v  ^  ^  and  on^y  if  ♦(t)  ■  exp(-t/2) ,  i.e.  x  NQ(r.In). 

Further,  if  x  ^  ECn(0,E,<|>)  with  E  >  0  and  p(x«0)  ■  0,  and  x  and 
s2  are  defined  by  (6.15),  by  the  corollary  of  Theorem  1  s2  *\>  G2(^^-;  y;  <t) 
if  and  only  if  E  is  a  generalized  inverse  of  (I-M) .  Let  T  be  an  orthogonal 
matrix  with  the  first  row  (1/v^T, . . .  tl/r/n) ,  then 

'0  0 

r'(i-M)r  - 

~  ~  ~  ~  0  1. 

-n-1 
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or 


i-m  »  r 


f  0  0  ] 


9  In- 1 


It  can  be  verified  that  any  generalized  inverse  (I-M)  of  (I-M)  can  be 
written  as 


(i-m)  -  r 


C11 

£12 

£21 

~n-lj 

rf  , 


where  c^f  c^  and  C21  are  arbitrary*  Hence  Z  must  have  the  above 
structure,  but  here  Z  >  0,  (c^)  Have  to  satisfy  the  conditions: 


and  C11  >  £12321* 


Example  2.  Assume  x  ^  EC^ ( 0 , , 4>)  ,  p(x»0)  -  0,  and  x  and  s  are 
defined  by  (6.15),  we  want  to  point  out  that 


(6.16) 


t  -  v'n(n-l)  \  *  t(n_1}  . 


Let  T  be  the  same  as  Example  1  and  /  ■  Tx,  then  y  ECft(0,In,<)>)  and 
y^  ■  r/n  x,  s^  -  y^.  If  z  ^  Nn(0,In),  then 


t  - 


/n-1  y1 
fn  2 

h  yi 


»  /n-l 


H'hW  i  ^  _  rr  _*X 


JTri'hf  JF% iiiii5 


- 


fn  2 

f|  Zi 


'v  t 


(n-l)  ’ 


because  (y^Ryll . yn/Ijj||)  *  (*1/|l*ll . hence  their  Borel  fu“ction® 

have  the  same  distribution. 
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Example  3.  Multivariate  t-distribution.  If  x  Nn(y,Z),  a  % 

2  2 

(i.e.,  8  ^  xv)  and  x  is  independent  of  8,  the  distribution  of 
y  =  A7  x/s  is  called  the  multivariate  t-distrlbutlon.  (cf.  Johnson  and 
Kotz  (1972)) .  When  n  ■  0,  E  ■  1^,  the  density  of  y  is 


(6.17) 


py<y> 


r(j  (n-K»)) 

(irv)^“r(i  v) 


(1  +  v 


-1 


y2j-5s(n+v) 


It  can  be  found  that  the  density  of  |  y  is 


(6.18) 


?lly| 


(u) 


2  rq  (n+v)) 


r(|  n)r(|  v)v 


n~l/i  ,  ~ 1  2. -*S(n+V)  . 

1  ..s^  U  (1  +  V  u  }  ,  u  >  0  . 


Obviously  the  multivariate  t-distribution  is  an  elliptically  contoured  dis¬ 
tribution.  If  x  'v-  Nn(p,E),  from  (1.2) 

x  $  y  +  RA’u(n) 

where  R  **»  ^  Is  independent  of  u^n^  and  A’ A  »  E.  Thus  if  y  «  0 


(6.19)  y  *  (AT  “)A'u^  , 

which  is  the  stochastic  representation  of  y.  The  density  (6.18)  of  /v  R/s 

'V  2  2 

can  be  motivated  by  the  F-distribution  because  —  (/v  ~)  ■  —  —  *\>  F(n,v). 

n  s  n  l 
s 

If  A  is  an  n  x  n  symmetric  matrix, 

x’Ax 

y’Ay  -  v  . 
s 

2  2 

As  x’Ax  is  independent  of  s  and  x’Ax  ^  If  and  only  if  AZA  -  A  and 

•  rk(A)  ■  k,  hence  r-  y’Ay  ^  F(k,v)  if  and  only  if  AZA  ■  A  and  rk(A)  ■  k. 
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Further  if  E  ■  IQ,  A  and  B  are  projection  matricea  with  rk(A)  ■  k, 
rk(B)  ■  l,  then  ^  'Ay,  y  'By)  is  distributed  according  to  multivariate 
inverted  Dirichlet  distribution  (cf.  p.  238,  Johnson  and  Kotz  (1972)). 

When  v  ■  1  the  multivariate  t -distribution  reduces  to  multivariate 
Cauchy  distribution. 


Example  A.  Uniform  distribution  in  the  unit  sphere  in  ]Rn.  If 
x  is  distributed  according  to  the  uniform  distribution  in  the  unit  sphere 
in  ]Rn,  its  density  is 


(6.20) 


0  otherwise  . 


It  can  be  shown  that  the  density  of  jjxjj  is 


P,|x,,(u>  •  nun“1  for  0  <  u  <  1  , 


and  the  density  of  |[x|p  is 


p  „(u)  *  for  0  <  u  <  1  . 

IWI2  2  “  " 


In  thin  cam  th.  danaity  of  Cjtoi  n~^ ;  il  la  (cf.  (4.1)) 
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“Vl#  x2'~  r  r-(n-2>(r2-*) 


n-1 

r  dr 


n  r(2)  2  ”  1  1  .  2  v  2  1 

It  n— k  *  n— 1  ^  "*) 

rcj)!*^)  n  1  ^ 


r(~) 


k_  ,  n+2-k 

x2  (1-x)  2 


0  <  x  <  1  . 


k.  n+2-k 

It  is  the  density  of  Bfc,  — = — ).  Similarly,  G  (a, , . . .  ,a  ,  ;a  ;<(»)  is 

l  l  mi  m— l  m 

equal  to  . ain-l'amfl^ *  where  2(a^  +  ••-.+ <»m)  ■  n. 

By  the  theorems  in  Section  4  we  obtain  the  following  conclusions:  If  A 

If  «yX9wk 

is  an  n  x  n  symmetric  matrix,  then  x'Ax  ^  BOy,  — j — )  if  and  only  if 
2 

A  “  A  and  rk(A)  »  k.  If  A^,  i  ”  l,...,m-l,  are  n  x  n  symmetric 


matrices,  then 


«  u  i  «■ — sa. 

(x'A.,x . x'A  .x)  'k  D  . -f— »  - - 

**  — 1~  —  -m-i-  m  l  l 


n+2-<x,  -  •  •  •  —o 


if  and  only  if  A^  -  A^  rkCA^  -  o^,  i  -  l,...,m-l  and  A^A^  -  0  for  i  4  j 
Clearly,  we  can  extend  the  above  conclusions  to  the  uniform  distribution 


in  ellipsoid 
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